c/3 



Idempotent and tropical mathematics. 
Complexity of algorithms and interval analysis 

G. L. Litvinov 1 

Institute for Information Transmission Problems of Russian Acad. Sci. 
Nagornaya, 27-4-72, Moscow 117186 Russia, tel. 7-499-1273804 



Abstract 

A very brief introduction to tropical and idempotent mathematics is presented. 
Tropical mathematics can be treated as a result of a dequantization of the tra- 
ditional mathematics as the Planck constant tends to zero taking imaginary 
values. In the framework of idempotent mathematics usually constructions and 
algorithms are more simple with respect to their traditional analogs. We espe- 
cially examine algorithms of tropical/idempotent mathematics generated by a 
collection of basic semiring (or semiheld) operations and other "good" opera- 
tions. Every algorithm of this type has an interval version. The complexity of 
this interval version coincides with the complexity of the initial algorithm. The 
interval version of an algorithm of this type gives exact interval estimates for the 
corresponding output data. Algorithms of linear algebra over idempotent and 
semirings are examined. In this case, basic algorithms are polynomial as well as 
their interval versions. This situation is very different from the traditional linear 
algebra, where basic algorithms are polynomial but the corresponding interval 
versions are NP-hard and interval estimates are not exact. 
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1. Introduction 

Tropical mathematics can be treated as a result of a dequantization of the 
traditional mathematics as the Planck constant tends to zero taking imaginary 
values. This kind of dequantization is known as the Maslov dequantization and 
it leads to a mathematics over tropical algebras like the max-plus algebra. The 
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Figure 1: Relations between idempotent and traditional mathematics. 



so-called idempotent dequantization is a generalization of the Maslov dequanti- 
zation. The idempotent dequantization leads to mathematics over idempotent 
semirings (exact definitions see below in sections 2 and 3). For example, the 
field of real or complex numbers can be treated as a quantum object whereas 
idempotent semirings can be examined as "classical" or "semiclassical" objects 
(a semiring is called idempotent if the semiring addition is idempotent, i.e. 
x®x = x), see [2H lil [27l [28l I2H . 

Tropical algebras are idempotent semirings (and semifields). Thus tropical 
mathematics is a part of idempotent mathematics. Tropical algebraic geometry 
can be treated as a result of the Maslov dequantization a ppli ed to the traditional 



algebraic geometry (O. Viro, G. Mikhalkin), see, e.g., [2ll. l52t l53l. l66l. 1671 168|. 
There are interesting relations and applications to the traditional convex geom- 
etry. 

In the spirit of N.Bohr's correspondence principle there is a (heuristic) cor- 
respondence between important, useful, and interesting constructions and re- 
sults over fields and similar results over idempotent semirings. A systematic 
application of this correspondence principle (which is a basic paradigm in idem- 
potcnt /tropical mathematcs) leads to a variety of theoretical and applied re- 
sults [II, M, 0, 0, HI EJ III , see Fig. 1 . 

The history of the subject is discussed, e.g., in [25(. There is a large list of 
references. 

In the framework of idempotent mathematics usually constructions and al- 
gorithms are more simple with respect to their traditional analogs (however, 
there exist NP-hard problems in tropical linear algebra). We especially exam- 
ine algorithms of tropical/idcmpotent mathematics generated by a collection of 
basic semiring (or semifield) operations and other "good" operations. Every al- 
gorithm of this type has an interval version. The complexity of this interval ver- 
sion coincides with the complexity of the initial algorithm. The interval version 
of an algorithm of this type gives exact interval estimates for the corresponding 
output data. Algorithms of linear algebra over idempotent and semirings are 
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examined. In this case, basic algorithms are polynomial as well as their interval 
versions. This situation is very different from the traditional linear algebra, 
where basic algorithms are polynomial but the corresponding interval versions 
are NP-hard and interval estimates are not exact. 

2. The Maslov dequantization 

Let R and C be the fields of real and complex numbers. The so-called max- 
plus algebra R max = RU {—00} is defined by the operations x © y = max{x, y] 
and x y = x + y. 

The max-plus algebra can be treated as a result of the Maslov dequantization 
of the semificld R+ of all nonnegative numbers with the usual arithmetics. The 
change of variables 

x n> u — h log x, 

where h > 0, defines a map $^ : R+ ->RU {—00}. This logarithmic transform 
was used by many authors. Let the addition and multiplication operations be 
mapped from R + to R U { — oc} by $>h, i-e. let 

u ®n v — ft.log(exp(u/ft.) + exp(u//i)), u v = u + v, 
= -00 = $ h (0), 1 = = $,,(!)• 

It can easily be checked that uffii, v — > max{w, v} as h — ► 0. Thus we get the 
semificld R ma x (i.e. the max-plus algebra) with zero = — 00 and unit 1 = as 
a result of this deformation of the algebraic structure in R+ . 

The semificld R m ax is a typical example of an idempotent semiring; this is 
a semiring with idempotent addition, i.e. x (33 x — x for arbitrary element x of 
this semiring. 

The semifield R max is also called a tropical algebra.The semifield R^J = 
<I > ; l (R + ) with operations (Bh and (i.e.+) is called a subtropical algebra. 

The semificld R m i n = R U {+00} with operations © = min and = + 
(0 = +00, 1 = 0) is isomorphic to R max - 

The analogy with quantization is obvious; the parameter h plays the role of 
the Planck constant. The map x <— > \x\ and the Maslov dequantization for R + 
give us a natural transition from the field C (or R) to the max-plus algebra 
Rmax- We will also call this transition the Maslov dequantization. In fact the 
Maslov dequantization corresponds to the usual Schrodinger dequantization but 
for imaginary values of the Planck constant (see below). The transition from 
numerical fields to the max-plus algebra R max (or similar semifields) in mathe- 
matical constructions and results generates the so called tropical mathematics. 
The so-called idempotent dequantization is a generalization of the Maslov de- 
quantization; this is the transition from basic fields to idempotent semirings in 
mathematical constructions and results without any deformation. The idempo- 
tent dequantization generates the so-called idempotent mathematics, i.e. math- 
ematics over idempotent semifields and semirings. Recently new versions of the 
Maslov dequantization appeared, see, e.g. (69l |. 
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Remark. The term 'tropical' appeared in [54| for a discrete version of the 
max-plus algebra (as a suggestion of Christian Choffrut). On the other hand 
V.P. Maslov used this term in 80s in his talks and works on economical applica- 
tions of his idempotent analysis (related to colonial politics). For the most part 
of modern authors, 'tropical' means 'over R max (or R m i n )' and tropical algebras 
are R max and R m i n . The terms 'max-plus', 'max-algebra' and 'min-plus' are 
often used in the same sense. 



3. Semirings and semifields 

Consider a set S equipped with two algebraic operations: addition © and 
multiplication 0. It is a semiring if the following conditions are satisfied: 

• the addition © and the multiplication are associative; 

the addition is commutative; 

the multiplication is distributive with respect to the addition ©: 

xQ(y®z) = (xQy)®(xQz) 

and 

(x © y) z = (x z) © (y z) 

for all x,y,z G S. 

A unity of a semiring S is an element 1 e S such that 10i = i01 = i for 
all x £ S. A zero of a semiring S is an element (if it exists) G S such that 
0/1 and © x = x, x = x = for all x £ S. A semiring S is called 
an idempotent semiring if x © x = x for all x £ S. A semiring S with a neutral 
element 1 is called a semifield if every nonzero element of S is invertible with 
respect to the multiplication. The theory of semirings and semifields is treated, 



e.g., in (14 1 . 



4. Idempotent analysis 

Idempotent analysis deals with functions taking their values in an idempo- 
tent semiring and the corresponding function spaces. Idempotent analysis was 
initially constructed by V. P. Maslov and his collaborators and then developed 
by many authors. The subject is presented in the book of V. N. Kolokoltsov 
and V. P. Maslov [22} (a version of this book in Russian was published in 1994). 

Let S be an arbitrary semiring with idempotent addition © (which is always 
assumed to be commutative), multiplication 0, and unit 1. The set S is supplied 
with the standard partial order ^<: by definition, a ^ b if and only if a © b = b. 
If the zero element exists, then all elements of S are nonnegative: ^ a for all 
a € S. Due to the existence of this order, idempotent analysis is closely related 
to the lattice theory, theory of vector lattices, and theory of ordered spaces. 
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Moreover, this partial order allows to model a number of basic "topological" 
concepts and results of idempotent analysis at t he p urel y algebraic level; this 



line of reasoning was examined systematically in [25|- [44| and [1CJ • 

Calculus deals mainly with functions whose values are numbers. The idem- 
potent analog of a numerical function is a map X — > S, where X is an arbitrary 
set and S is an idempotent semiring. Functions with values in S can be added, 
multiplied by each other, and multiplied by elements of S pointwise. 

The idempotent analog of a linear functional space is a set of S- valued func- 
tions that is closed under addition of functions and multiplication of functions by 
elements of S, or an S-semimodule. Consider, e.g., the S'-semimodule B(X, S) 
of all functions X — > S that are bounded in the sense of the standard order on 
S. 

If S = R-max, then the idempotent analog of integration is defined by the 
formula 



Hp) — / <P( X ) dx — sup f(x), (1) 
Jx xex 

where f £ B(X, S). Indeed, a Riemann sum of the form ^2 f{xi)-<Ji corresponds 

i 

to the expression @ y(xj)©<7j = mait{if(xi) + <Ji}, which tends to the right-hand 

i ' 

side of (1) as <Xj — > 0. Of course, this is a purely heuristic argument. 

Formula (1) defines the idempotent (or Maslov) integral not only for func- 
tions taking values in R ma x, but also in the general case when any of bounded 
(from above) subsets of S has the least upper bound. 

An idempotent (or Maslov) measure on X is defined by the formula m^,(Y) = 
supip(x), where tp € B(X,S) is a fixed function. The integral with respect to 

this measure is defined by the formula 



Iipif) = / f{x)dm^,= I f(x)Giilj(x)dx=sup(f(x)Qip(x)). (2) 
Jx Jx xex 

Obviously, if S = R m in, then the standard order is opposite to the conven- 
tional order <, so in this case equation (2) assumes the form 

f{x)dm^= I f(x)Qip(x)dx= inf (f(x) ip(x)), 
x Jx xeX 

where inf is understood in the sense of the conventional order < . 



5. The superposition principle and linear problems 

Basic equations of quantum theory are linear; this is the superposition princi- 
ple in quantum mechanics. The Hamilton-Jacobi equation, the basic equation 
of classical mechanics, is nonlinear in the conventional sense. However, it is 
linear over the semirings R max and R m i n - Similarly, different versions of the 
Bellman equation, the basic equation of optimization theory, are linear over 
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suitable idempotent semirings; this is V. P. Maslov's idempotent superposition 
principle, see 4|| 49|, 5(j. For instance, the finite-dimensional stationary Bell- 



man equation can be written in the form X = H © X © F, where X, H, F 
are matrices with coefficients in an idempotent semiring S and the unknown 
matrix X is determined by H and F [7, 8]. In particular, standard problems of 
dynamic programming and the well-known shortest path problem correspond 
to the cases S = R max and S = R m in, respectively. It is known that principal 
optimization algorithms for finite graphs correspond to standard methods for 
solving systems of linear equations of this type (i.e., over semirings). Specif- 
ically, Bellman's shortest path algorithm corresponds to a version of Jacobi's 
algorithm, Ford's algorithm corresponds to the Gauss-Seidel iterative scheme, 
etc. 0,13. 

The linearity of the Hamilton- Jacobi equation over R m j n and R max , which 
is the result of the Maslov dequantization of the Schrodinger equation, is closely 
related to the (conventional) linearity of the Schrodinger equation and can be 
deduced from this linearity. Thus, it is possible to borrow standard ideas and 
methods of linear analysis and apply them to a new area. 

Consider a classical dynamical system specified by the Hamiltonian 

JV 2 
i—l 1 

where x — {x\, . . . , xn) are generalized coordinates, p — (p%, . . . ,pn) are gener- 
alized momenta, m, are generalized masses, and V(x) is the potential. In this 
case the Lagrangian L(x, x, t) has the form 

N .0 



L(x,x,t) = J2 mi ^T ~ V W 



i=l 

where x — (ii, . . . , xjsr), Xi = dxi/dt. The value function S(x, t) of the action 
functional has the form 

S= [ L(x(t),x(t),t)dt, 
J t 

where the integration is performed along the factual trajectory of the system. 
The classical equations of motion are derived as the stationarity conditions for 
the action functional (the Hamilton principle, or the least action principle). 

For fixed values of t and to and arbitrary trajectories x(t), the action func- 
tional S = S(x(t)) can be considered as a function taking the set of curves 
(trajectories) to the set of real numbers which can be treated as elements of 
R m i n . In this case the minimum of the action functional can be viewed as the 
Maslov integral of this function over the set of trajectories or an idempotent 
analog of the Euclidean version of the Feynman path integral. The minimum 
of the action functional corresponds to the maximum of e~ s , i.e. idempotent 
integral J^ aths j e~ s ^ x ^ D{x(t)} with respect to the max-plus algebra R m ax- 
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Thus the least action principle can be considered as an idcmpotent version of 
the well-known Feynman approach to quantum mechanics. The representation 
of a solution to the Schrodinger equation in terms of the Feynman integral cor- 
responds to the Lax-Olemik solution formula for the Hamilton-Jacobi equation. 

Since dS/dxi = Pi, dS/dt = —H(p, x), the following Hamilton- Jacobi equa- 
tion holds: 

dS „ ( OS 



Quantization leads to the Schrodinger equation 

~^ = H 1 /, = H(p i ,x i )^, (4) 

where ip = tp(x,t) is the wave function, i.e., a time-dependent element of the 
Hilbert space L 2 (R N ), and H is the energy operator obtained by substitution of 
the momentum operators Pi = jg§- and the coordinate operators x~i : tp i-> Xiip 
for the variables pi and Xi in the Hamiltonian function, respectively. This equa- 
tion is linear in the conventional sense (the quantum superposition principle). 
The standard procedure of limit transition from the Schrodinger equation to the 
Hamilton-Jacobi equation is to use the following ansatz for the wave function: 
ip(x,t) = a(x,t)e lS ( x ' t ^ ri , and to keep only the leading order as h — > (the 
'semiclassical' limit). 

Instead of doing this, we switch to imaginary values of the Planck constant h 
by the substitution h = ih, assuming h > 0. Thus the Schrodinger equation (4) 
turns to an analog of the heat equation: 

h ^7 = H (- h 4r^i) u ^ ( 5 ) 



dt \ dx. 

where the real-valued function u corresponds to the wave function ip. A similar 
idea (the switch to imaginary time) is used in the Euclidean quantum field 
theory; let us remember that time and energy are dual quantities. 

Linearity of equation (4) implies linearity of equation (5). Thus if u\ and ui 
are solutions of (5), then so is their linear combination 

M = A 1 M 1 +A 2 M 2 - (6) 

Let S = h\nu or u — e s l h as in Section 2 above. It can easily be checked 
that equation (5) thus turns to 

dS . . A 1 fdS\ 2 1 8 2 S 

Thus we have a transition from (4) to (7) by means of the change of variables 
ip = e s/h . Note that |V>| = e ReS/h , where RcS is the real part of S. Now let us 
consider S as a real variable. The equation (7) is nonlinear in the conventional 
sense. However, if Si and ^2 are its solutions, then so is the function 
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obtained from (6) by means of our substitution S — Mnu. Here the general- 
ized multiplication coincides with the ordinary addition and the generalized 
addition (Bh is the image of the conventional addition under the above change 
of variables. As h — > 0, we obtain the operations of the idempotent semiring 
R max , i.e., © = max and = +, and equation (7) turns to the Hamilton-Jacobi 
equation (3), since the third term in the right-hand side of equation (7) vanishes. 

Thus it is natural to consider the limit function S = Ai Si © A2 S2 as 
a solution of the Hamilton-Jacobi equation and to expect that this equation 
can be treated as linear over R max - This argument (clearly, a heuristic one) 
can be extended to equations of a more general form. For a rigorous treatment 
of (semiring) linearity for these equations see, e.g., [22|,[30j]. Notice that if h is 
changed to —h, then we have that the resulting Hamilton-Jacobi equation is 
linear over R m i n . 

The idempotent superposition principle indicates that there exist important 
nonlinear (in the traditional sense) problems that are linear over idempotent 
semirings. The idempotent linear functional analysis (see [H Hi m, 
35l 3rl 43l H . 51 , H 16 1 ) is a natural tool for investigation of those nonlinear 
infinite-dimensional problems that possess this property. In practice infinite- 
dimensional problems can be approximated by finite-dimensional problems. So 
algorithms of idempotent linear algebras are especially important (because of 
the superposition principle!). Below some examples are examined. 



6. Applications 

There are very many important applications of tropical/idempotent math- 
ematics (and especially the correspondence and superposition principles) in- 
cluding optimization and control, algebraic geometry, dynamic programming, 
differential equations, mathematical biology, mathematical physics and chem- 
istry, transport and energoenergetic netwoks, interval analysis, mathematical 
economics, game theory, computer technology etc., see, e.g. 

Bum mm 

[13, HI, H M, 0, IH EE 0, SS EH HI Hi Hi HI, |^| . Applications of 
the idempotent correspondence principles to software and hardware desig n are 
examined, e.g. in 27, 28, 38|. Some applications are discucced in [32I [3lf. 



7. Positive semirings and basic operations 

7.1. Some definitions 

Let the semiring S be partially ordered (see, e.g. [1, HH and Subsection 9.2 
below) by a relation ^< such that is the least element and the inequality x -< y 
implies that x®z<y®z,xQz<yQz, and z x ^ z y for all x,y, z G S; 
in this case the semiring S is called positive (see, e.g., [Hj]). 

Recall that a semiring S is called idempotent if x © x = x for all x G S. 
In this case the addition © defines a canonical (or standard) partial order ^ 
on the semiring S by the rule: x < y iff x © y — y. It is easy to prove that 
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any idcmpotcnt semiring is positive with respect to this order. Note also that 
x © y — sup{x, y} with respect to the canonical order. In the sequel, we shall 
assume that all idempotent semirings are ordered by the canonical partial order 
relation. 

We shall say that a positive (e.g., idempotent) semiring S is complete if it is 
complete as an ordered set. This means that for every subset T C S there exist 
elements sup T £ S and inf T £ S. 

The most well-known and important examples of positive semirings are "nu- 
merical" semirings consisting of (a subset of) real numbers and ordered by the 
usual linear order < on R: the semiring R + with the usual operations © = +, 
= • and neutral elements = 0, 1 = 1, the semiring R max = R U {— oo} 
with the operations © = max, = + and neutral elements = — oo, 1=0, 
the semiring R max = Rmax U {oo}, where x ^ oo, x © oo = oo for all x, 
xGog = oo0a; = oo if x ^ 0, and O0oo = oo0O, the semirings Z max = ZU{— oo} 
and Z max = Z max U {+00} (subsemirings in R max and R max ), and the semir- 
ing <S^ ax min = [a, b], where — oo<a<6<+oo, with the operations © = max, 
= min and neutral elements = a, 1 = b. The semirings R maX5 Rmax, 
Z max , Z max and S^ a ^ min = [a, 6] are idempotent. The semirings R ma x, Z max , 

■Smax min , R + = R + U{°°} are complete. Remind that every partially ordered 
set can be imbedded to its completion (a minimal complete set containing the 
initial one). We shall say that all these semirings (as well as algebras isomorphic 
to them) are basic numerical semiring. These semirings are complete or their 
completions are complete semirings. 

Denote by WL&t mn (S) a set of all matrices A = (a^) with m rows and 
n columns whose coefficients belong to a semiring S. The sum A © B of 
matrices A,B £ Mat mn (5 l ) and the product AB of matrices A £ Matj m (if>) 
and B £ Mat mn (S) are defined according to the usual rules of linear algebra: 
A © B = (a,ij © bij) £ Mat mn (5) and 



where A £ Matj m (S) and B £ M&t mn (S). Note that we write AB instead of 
A 5. 

If the semiring S is positive, then the set Mat mn (5) is ordered by the relation 
A = (aij) ^ B = (bij) iff ^ bij in S for all 1 < i < m, 1 < j < n. 

The matrix multiplication is consistent with the order ^ in the following 
sense: if A, A' £ Mat im (S), B,B' £ Mat„ m (5) and A ^ A', B ^ B', then 
AB ^ A'B' in Mat;„(S l ). The set Mat nn (S) of square (n x n) matrices over a 
[positive, idempotent] semiring S forms a [positive, idempotent] semiring with 
a zero element O = (o^), where = 0, 1 < i,j < n, and a unit element 
J = (<%), where Sij = 1 if i = j and 5ij — otherwise. 

The set Mat„ n is an example of a noncommutative semiring if n > 1. 
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7.2. Closure operations 

Let a positive semiring S be endowed with a partial unary closure operation 
or Kleene star operation * such that a ^ b implies a* ^ b* and a* — l©(a*0a) = 
1 © (a a*) on its domain of definition. In particular, 0* = 1 by definition. 

These axioms imply that a* = 1 © a © a 2 © ■ • • © (a* a n ) if n > 1. Thus 
x* can be considered as a 'regularized sum' of the series a* = 1 © a © a 2 © ... . 
In a positive semiring, provided that it is closed under taking bounded ordered 
sup-operations and the operations © and distribute over such sup-operations 
we can define 

a* :=supl©affi...©a fc , (7.1) 

fc>0 

if the sequence on the r.h.s. is bounded. In this case a* is the least solution 

of the equations x = ax © 1 and x = xa © 1, and more generally a*b is the the 
least solution of the Bellman equations x = ax © b and x = xa © b. So if S is 
complete, then the closure operation is well-defined for every element x 6 S. 
In the case of idempotcnt addition (|7.1[) becomes particularly nice: 

a*=0a , =supa ! . (7.2) 

1 ii numerical semirings the operation * is usually very easy to implement: 
x* = (1 — x) if x < 1 in R+, or R + and x* = oo if x > 1 in R + ; x* = 1 
if x ^ 1 in R m ax and R max , x* = cx) if x ^ 1 in R ma x, x* = 1 for all x in 
Sl?' a b } In all other cases x* is undefined. 

The closure operation in matrix semirings over a positive semiring S can be 
defined inductively: A* = (an)* = (a| x ) in Matn(S') and for any integer n > 1 
and any matrix 

'An A 12 ' 
A 2 \ A12. 



A 



where An e Mat fefe (S), A 12 £ Mat fe „_ fc (S'), A 2 i S Mat„_ fcfc (5), A 22 £ Mat„_ feri _ 
1 < fc < n, by defintion, 

A4*i © A* lx A 12 D*A 2 xA* n A* n A 12 D*\ 
A*=[ , (7.3) 

V ^21^*! D* J 

where D = A 22 © A 2 iAl 1 Ai2- It can be proved that this definition of A* implies 
that the equalities A* = A* A © / = A A* © 7" are satisfied, and thus A* is a 
'regularized sum' of the series I © A © A 2 © ... . Moreover, in the case when A* 
is defined as the least solution of A* = A* A © 7 and A* = yL4* © I, it can be 
shown that it satisfies (|7.3j) . 

Note that this recurrence relation coincides with the formulas of escalator 
method of matrix inversion in the traditional linear algebra over the field of real 
or complex numbers, up to the algebraic operations used. Hence this algorithm 
of matrix closure requires a polynomial number of operations in n, see 0, [Til 
M EH, S, 0] for more details. 
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Let S be a positive semiring. The discrete stationary Bellman equation or 
matrix Bellman equation has the form 

X = AX®B, (7.4) 

where A € MeA, nn (S), X, B £ Mat ns (5), and the matrix X is unknown. Let A* 
be the closure of the matrix A. It follows from the identity A* = A* A © / that 
the matrix A*B satisfies this equation. As in the scalar case, it can be shown 
that for positive semirings under reasonable distributivity assumptions, if A* is 
defined as in (|7.1I) then A*B is the least in the set of solutions to equation (|7.4p 
with respect to the partial order in Mat„ s (5). Recall that in the idempotent 
case 

A* = @A i = supA\ (7.5) 

*>o ^° 

Consider also the case when A = (ay) is nxn strictly upper-triangular (such 
that ay = for i > j), or n X n strictly lower-triangular (such that a™ = 
for i < j). In this case A™ = O, the all-zeros matrix, and it can be shown by 
iterating X = AX © / that this equation has a unique solution, namely 

A* = I®A®...®A n - 1 . (7.6) 

Curiously enough, formula (|7.6j) works more generally in the case of numerical 
idempotent semiring R max (and other idempotent semirings): in fact, the se- 
ries (|7.5|) converges there if and only if it can be truncated to (|7.6|) . This is 
closely related to the principal path interpretation of A* explained in the next 
subsection. 

7.3. Weighted directed graphs and matrices over semirings 

Suppose that S is a semiring with zero and unity 1. It is well-known that 
any square matrix A = (a^ ) E Mat nn (S l ) specifies a weighted directed graph. 
This geometrical construction includes three kinds of objects: the set X of n 
elements x\,...,x n called nodes, the set T of all ordered pairs (xi,Xj) such that 
ciij 7^ called arcs, and the mapping A: T — ¥ S such that A(xi,Xj) — aij. The 
elements dy of the semiring S are called weights of the arcs. 

Conversely, any given weighted directed graph with n nodes specifies a 
unique matrix A S Mat nra (S l ). 

This definition allows for some pairs of nodes to be disconnected if the corre- 
sponding element of the matrix A is and for some channels to be "loops" with 
coincident ends if the matrix A has nonzero diagonal elements. This concept is 
convenient for analysis of parallel and distributed computations and design of 
computing media and networks. 

Recall that a sequence of nodes of the form 

P = (yo,2/i, •••,2/fe) 

with k > and (j/j, 2/i+i) £ L, i — 0, .. . ,k — 1, is called a path of length k 
connecting yo with yu- Denote the set of all such paths by Pk(yo,yk)- The 
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weight A(p) of a path p S P k {yo, Uk) is defined to be the product of weights of 
arcs connecting consecutive nodes of the path: 



[fc] 



A{ P ) = A(y , yi) • • • © A{y k - X ,y k ). 

By definition, for a 'path' p G Po(xi, Xj) of length k — the weight is 1 if i = j 
and otherwise. 

For each matrix A s Mat nrl (5 l ) define A° = I = (Sij) (where Sij = 1 if i = j 
and Sij — otherwise) and A k = AA k ^ 1 , k > 1. Let be the (i,j)th element 
of the matrix A k . It is easily checked that 

^3 S'l • • • a i k -!i k ■ 

io=i, ik=j 
l<ii,...,ifc_i <n 

Thus afj is the supremum of the set of weights corresponding to all paths of 
length k connecting the node Xi = Xi with Xi k — Xj . 

Let A* be defined as in (I7.5[) . Denote the elements of the matrix A* by a*j , 
i,j = l,...,n; then 

4= 

0<fc<oo p£P k (xi,Xj) 

The closure matrix A* solves the well-known algebraic path problem, which 
is formulated as follows: for each pair (:rj, xj) calculate the supremum of weights 
of all paths (of arbitrary length) connecting node Xi with node Xj . The closure 
operation in matrix semirings has been studied extensively (see, e.g., 0, H, O 



U, |15|, ll3 [22J, |46| and references therein). 



Example 7.1 (The shortest path problem). Let S — R m i n , so the weights are 
real numbers. In this case 

A(p) = A(y , yi) + A( yi ,y 2 ) + ■■■ + A{y k ^y k ). 

If the element specifies the length of the arc (xi,Xj) in some metric, then 
a*j is the length of the shortest path connecting Xi with Xj . 

Example 7.2 (The maximal path width problem). Let S = R U {0, 1} with 
® = max, = min. Then 



a,. = max 
J pe U Pkixuxj) 

k>l 



A(p), A(p) = mm(A(y , y x ), A(y k -i,y k )). 



If the element specifies the "width" of the j), then the width of a 

path p is defined as the minimal width of its constituting arcs and the element 
a*j gives the supremum of possible widths of all paths connecting xi with Xj . 

Example 7.3 (A simple dynamic programming problem). Let S = R max and 
suppose dij gives the profit corresponding to the transition from Xi to Xj . De- 
fine the vector B = (pi) €E Mat n i(R max ) whose element hi gives the terminal 
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profit corresponding to exiting from the graph through the node Xi. Of course, 
negative profits (or, rather, losses) are allowed. Let m be the total profit corre- 
sponding to a path p G Pp.{xi,Xj), i.e. 

m = A(p) + bj. 

Then it is easy to check that the supremum of profits that can be achieved on 
paths of length k beginning at the node Xi is equal to (A k B)i and the supremum 
of profits achievable without a restriction on the length of a path equals (A*B)i. 

Example 7.4 (The matrix inversion problem). Note that in the formulas of 
this section we are using distributivity of the multiplication with respect to 
the addition © but do not use the idempotency axiom. Thus the algebraic 
path problem can be posed for a nonidempotent semiring S as well (this is 
well-known). For instance, if S = R, then 

A* = I + A + A 2 + ■ ■ ■ = (I - A)' 1 . 

If 1 1 A 1 1 > 1 but the matrix I — A is invertible, then this expression defines a 
regularized sum of the divergent matrix power series J2i>o ^ ' ■ 

We emphasize that this connection between the matrix closure operation 
and solution to the Bellman equation gives rise to a number of different algo- 
rithms for numerical calculation of the closure matrix. All these algorithms are 
adaptations of the well-known algorithms of the traditional computational linear 
algebra, such as the Gauss- Jordan elimination, various iterative and escalator 
schemes, etc. This is a special case of the idempotent superposition principle. 

In fact, the theory of the discrete stationary Bellman equation can be de- 
veloped using the identity A* = AA* © I as an additional axiom without any 
substantive interpretation (the so-called closed semirings, see, e.g., 



Such theory can be based on the following identities, true both for the case 
of idempotent semirings with path interpretation, and the real numbers with 
conventional arithmetic (assumed that A and B have appropriate sizes): 



(A®B)* = (A*B)*A* 
(AB)*A — A(BA)* . 



(7.7) 



7. 4- Basic operations 

Let S be an idempotent or positive semiring. Then S is a partial ordered 
set (or poset) with respect to the canonical order. Suppose that S is a lattice, 
i.e. for each pair of elements x, y, there exists the least lower bound x V y called 
supremum and the greatest lower bound x A y called infimum. See details in 
Subsection 9.1 below. For basic numerical semirings (see Subsection 7.2 above) 
these operations are maximum and minimum. 

We shall say that the semiring operations, supremum, infinum and the unary 
closure operation (Kleene star-operation) are basic operations. 

We shall say that S is a completed semifield if S is a complete semiring and 
S without the element sup S is a semifield. Then the (unary) inversion x H> x^ 1 
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is obviously well defined for every element. If S is a semificld or completed 
semifield, then we shall say that the inversion operation is also basic. 

For basic numerical positive semirings all the basic operations are very simple 
and easy for computer implementations. 



8. Algorithms of tropical/idempotent mathematics 



8.1. The correspondence principle for computations 

Of course, the idempotent correspondence principle is valid for algorithms as 
well as for their software and hardware implementations [27, 
Thus: 



// we have an important and interesting numerical algorithm, then 
there is a good chance that its semiring analogs are important and 
interesting as well. 

In particular, according to the superposition principle, analogs of linear al- 
gebra algorithms are especially important. Note that numerical algorithms for 
standard infinite-dimensional linear problems over idempotent semirings (i.e., 
for problems related to idempotent integration, integral operators and transfor- 
mations, the Hamilton- Jacobi and generalized Bellman equations) deal with the 
corresponding finite-dimensional (or finite) "linear approximations" . Nonlinear 
algorithms often can be approximated by linear ones. Thus the idempotent 
linear algebra is a basis for the idempotent numerical analysis. 

Moreover, it is well-known that linear algebra algorithms easily lend them- 
selves to parallel computation; their idempotent analogs admit parallelization 
as well. Thus we obtain a systematic way of applying parallel computing to 
optimization problems. In this paper we do not deal with parallel algorithms 
and their implementations. 

8.2. Universal algorithms 

Computational algorithms are constructed on the basis of certain primitive 
operations. These operations manipulate data that describe "numbers." These 
"numbers" are elements of a "numerical domain," i.e., a mathematical object 
such as the field of real numbers, the ring of integers, numerical and idempotent 
semirings and semifields. 

In practice, elements of the numerical domains are replaced by their com- 
puter representations, i.e., by elements of certain finite models of these domains. 
Examples of models that can be conveniently used for computer representation 
of real numbers are provided by various modifications of floating point arith- 
metics, approximate arithmetics of rational numbers (4o| . and interval arith- 
metics. The difference between mathematical objects ("ideal" numbers) and 
their finite models (computer representations) results in computational (e.g., 
rounding) errors. 
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An algorithm is called universal if it is independent of a particular numer- 
ical domain and/or its computer representation. A typical example of a uni- 
versal algorithm is the computation of the scalar product (x, y) of two vectors 

x = (si, . . . , x n ) and y = . . . , y n ) by the formula (x, y) = x x y\ H V x n y n . 

This algorithm (formula) is independent of a particular domain and its com- 
puter implementation, since the formula is well-defined for any semiring. It is 
clear that one algorithm can be more universal than another. For example, the 
simplest Newton-Cotes formula, the rectangular rule, provides the most uni- 
versal algorithm for numerical integration. In particular, this formula is valid 
also for idempotent integration (that is, over any idempotent semiring, see e.g. 
22, IS IH 4i|[50j]. Other quadrature formulas (e.g., combined trapezoid rule 



or the Simpson formula) are independent of computer arithmetics and can be 
used (e.g., in the iterative form) for computations with arbitrary accuracy. In 
contrast, algorithms based on Gauss-Jacobi formulas are designed for fixed ac- 
curacy computations: they include constants (coefficients and nodes of these 
formulas) defined with fixed accuracy. (Certainly, algorithms of this type can 
be made more universal by including procedures for computing the constants; 
however, this results in an unjustified complication of the algorithms.) 

Modern achievements in software development and mathematics make us 
consider numerical algorithms and their classification from a new point of view. 
Conventional numerical algorithms are oriented to software (or hardware) im- 
plementation based on floating point arithmetic and fixed accuracy. However, 
it is often desirable to perform computations with variable (and arbitrary) ac- 
curacy. For this purpose, algorithms are required that are independent of the 
accuracy of computation and of the specific computer representation of numbers. 
In fact, many algorithms are independent not only of the computer represen- 
tation of numbers, but also of concrete mathematical (algebraic) operations on 
data. In this case, operations themselves may be considered as variables. Such 
algorithms are implemented in the form of generic programs based on abstract 
data types that are defined by the user in addition to the predefined types pro- 
vided by the language. The corresponding program tools appeared as early as in 
Simula-67, but modern object-oriented languages (like C ++ , see, e.g., [471 \6l\) 
are more convenient for generic programming. Computer algebra algorithms 
used in such systems as Mathematica, Maple, REDUCE, and others are also 
highly universal. 

A different form of universality is featured by iterative algorithms (begin- 
ning with the successive approximation method) for solving differential equa- 
tions (e.g., methods of Euler, Euler-Cauchy, Runge-Kutta, Adams, a number 
of important versions of the difference approximation method, and the like), 
methods for calculating elementary and some special functions based on the 
expansion in Taylor's series and continuous fractions (Pade approximations). 
These algorithms are independent of the computer representation of numbers. 

The concept of a generic program was introduced by many authors; for 
example, in 2J| such programs were called 'program schemes.' In this paper, 
we discuss universal algorithms implemented in the form of generic programs. 
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8.3. The correspondence principle for hardware design 

A systematic application of the correspondence principle to computer calcu- 
lations leads to a unifying approach to software and hardware design. 

The most important and standard numerical algorithms have many hardware 
realizations in the form of technical devices or special processors. These devices 
often can be used as prototypes for new hardware units generated by substitution 
of the usual arithmetic operations for its semiring analogs and by addition tools 
for performing neutral elements and 1 (the latter usually is not difficult). Of 
course, the case of numerical semirings consisting of real numbers (maybe except 
neutral elements) and semirings of numerical intervals is the most simple and 
natural . Note that for scmifields (including R max and R-min) the operation of 
division is also defined. 

Good and efficient technical ideas and decisions can be transposed from 
prototypes into new hardware units. Thus the correspondence principle gen- 
erated a regular heuristic method for hardware design. Note that to get a 
patent it is necessary to present the so-called 'invention formula', that is to 
indicate a p rot otyp e for the suggested device and the difference between these 
devices (23, M, El HI M, E| ■ 

Consider (as a typical example) the most popular and important algorithm 
of computing the scalar product of two vectors: 

(x, y) = xxyx + x 2 y 2 H hx„y n . (8.1) 

The universal version of (12) for any semiring A is obvious: 

(x, y) = (xi Vl ) © (x 2 y 2 ) © • • ■ © (x„ y n ). (8.2) 

In the case A = R max this formula turns into the following one: 

(x, y) = max{xi + yi, x 2 + 2/2, • • • ,x„ + y n }- (8.3) 

This calculation is standard for many optimization algorithms, so it is useful 
to construct a hardware unit for computing (14). There are many different 
devices (and patents) for computing (12) and every such device can be used as 
a prototype to construct a new device for computing (14) and even (13). Many 
processors for matrix multiplication and for other algorithms of linear algebra 
are based on computing scalar products and on the corresponding "elementary" 
devices respectively, etc. 

8.4- The correspondence principle for software design 

Software implementations for universal semiring algorithms are not as effi- 
cient as hardware ones (with respect to the computation speed) but they are 
much more flexible. Program modules can deal with abstract (and variable) 
operations and data types. Concrete values for these operations and data types 
can be defined by the corresponding input data. In this case concrete operations 
and data types are generated by means of additional program modules. For pro- 
grams written in this manner it is convenient to use special techniques of the 
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so-called object oriented (and functional) design, see, e.g., [47J,l6l|. Fortunately, 
powerful tools supporting the object-oriented software design have recently ap- 
peared including compilers for real and convenient programming languages (e.g. 
C ++ and Java) and modern computer algebra systems. 

Recently, this type of programming technique has been dubbed the so-called 
generic programming (see, e.g., [sH)- To help automate the generic program- 
ming, the so-called Standard Template Library (STL) was developed in the 



framework of C ++ [6l|, |64|. However, high-level tools, such as STL, possess 



both obvious advantages and some disadvantages and must be used with cau- 
tion. 



8.5. Complexity of algorithms in idempotent mathematics 

We shall use the well known standard terminology of the complexity the- 
ory (time complexity, space complexity, asymptotic computational complexity, 
polynomial complexity, NP-hard problems etc.), see, e.g., the corresponding 
Wikipedia articles. The time complexity of an algorithm quantifies the amount 
of time taken by this algorithm to run as a function of the size of the input to the 
problem. The time complexity of an algorithm is commonly expressed using big 
O notation, which suppresses multiplicative constants and lower order terms. 
When expressed this way, the time complexity is said to be described asymptot- 
ically, i.e., as the input size goes to infinity. The time complexity is commonly 
estimated by counting the number of elementary operations performed by the 
algorithm. In idempotent mathematics ( and mathematics over positive semir- 
ings) the elementary operations are basic operations described in Subsection 7.4 
above. 

For the space complexty, the situation is quite similar. The following def- 
inition is important for us. We shall say that two algorithms have the same 
complexity if they have the same asymptotic time and space complexity. 

In principle idempotent mathematcs and its algorithms are more simple 
with respect to traditional mathematics. That is why the most important al- 
gorithms of idempotent mathematics (and espcially idempotent linear algebra) 
are polynomial. For example, many algorithms of solving the stationary discrete 
(matrix) Bellman equations have the complexity 0(n 3 ), see 0, S, EH, H, S, [3^, 

olynomial algorithms of linear algebra are examined 
32, 65| . However, NP-hard problems exist (e.g., in 




tropical and idempotent linear algebra), see [l7, 18, 32, 62 1. Some frontiers of 
polynomial computations in tropical geometry are investigated in (65l |. 



9. Interval analysis in idempotent mathematics and complexity of 
algorithms 

9.1. Interval extensions of algorithms 

Interval analysis appears for treating input and output data with uncertain- 
ties (interval data). Traditional interval analysis is a nontrivial and popular 

mathematical area, see, e.g., p, m, m, 0, hi. An "idempotent" version of 
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interval analysis (and moreover interval analysis over positive semirings) ap- 
peared in 45 L 46 L 63 L Rather many publications on the subject appeared later, 
see, e.g., laElEalHSlHB. Interval analysis over the positive semiring R + was 



discussed in [3]. In the framework of idempotent mathematics, interval analysis 
gives exact interval solutions without any conditions of smallness on uncertainty 
intervals. 

Let a set S be partially ordered by a relation ^. Below (in Subsection 9.2 
partially ordered sets (or posets for the sake of brevity) will be discussed in 
details. A closed interval in S is a subset of the form x = [x, x] = { x E S | 
x ^ x ^ x}, where the elements x < x are called lower and upper bounds of 
the interval x. The order < induces a partial ordering on the set of all closed 
intervals in S: x < y iff x ^ y and x ^ y. 

A weak interval extension I(S) of a positive semiring S is the set of all 
closed intervals in S endowed with operations © and defined as x © y = 
[x © y , x © y] , xOy = [x0y,x0y] and a partial order induced by the order 
in S. The closure operation in I(S) is defined by x* = [x*,x*]. There are some 
other interval extensions (including the so-called strong interval extension |46j |) 
but the weak extension is more convenient. 

The extension I(S) is a positive semiring; I(S) is idempotent if S is an 
idempotent semiring. A universal algorithm over S can be applied to I(S) 
and we shall get an interval version of the initial algorithm. However, there 
are some conditions for interval extensions of algorithms to be sure that these 
extensions are good enough, see below. Usually both the versions have the same 
complexity. For the discrete stationary Bellman equation and the corresponding 
optimization problems on graphs, interval analysis was examined in 4a, III in 
details. Other problems of idempotent linear algebra were examined in 9j, Il3l 
II [HI, [57 1 . 

Idempotent mathematics appears to be remarkably simpler than its tradi- 
tional analog. For example, in traditional interval arithmetic, multiplication 
of intervals is not distributive with respect to addition of intervals, whereas 
in idempotent interval arithmetic this distributivity is preserved. Moreover, in 
traditional interval analysis the set of all square interval matrices of a given 
order does not form even a semigroup with respect to matrix multiplication: 
this operation is not associative since distributivity is lost in the traditional 
interval arithmetic. On the contrary, in the idempotent (and positive) case as- 
sociativity is preserved. Finally, in traditional interval analysis some problems 
of linear algebra, such as solution of a linear system of interval equations, can 
be very difficult (more precisely, the y a re iVP- hard, see [ll|,[l3, 23 [ and refer- 
ences therein). It was noticed in (4BI |46| that in the idempotent case solving an 
interval linear system requires a polynomial number of operations (similarly to 
the usual Gauss elimination algorithm). Two properties that make the idem- 
potent interval arithmetic so simple are monotonicity of arithmetic operations 
and positivity of all elements of an idempotent semiring. 

Usually interval estimates in idempotent mathematics are exact. In the 
traditional theory such estimates tend to be overly pessimistic. 
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9.2. Intervals in partially ordered sets and interval regular mappings 

Let us start with some basic notions of the theory of lattices and partially 
ordered sets. The reader is referred to Q for more information. 

Definition 9.1. Binary relation < on a set S is called a partial order if it 
satisfies the following axioms: 1) a < a, 2) a ■< b and b -< a imply a = b, 3) 
a <b and b < c imply a ■< c. In this case S is called a partially ordered set or, 
briefly, a poset. 

If S is a Cartesian product Si x S 2 where Si and S 2 are posets ordered 
with ^1 and < 2 respectively, one can naturally introduce relation ^ on S, by 
(xi,yi) r< (^2,2/2) xi ^1 x 2 and yi ^ 2 2/2- 

For the dynamics, consider mappings of partially ordered sets. A mapping 
<fi: S — > T is a morphism of partially ordered sets S and T if x ■< y implies 
4>{x) < 4>(y). This mapping is an isomorphism if it is a bijection (i.e. one-to-one 
correspondence) between S and T. Also note that if 4>\ is a morphism between 
Si and Ti, and 02 is a morphism (resp. isomorphism) between S 2 and T2, then 
mapping <f>i x (f> 2 : (x,jy) 1— > ((j)i(x),4) 2 (y)) is a morphism (resp. isomorphism) 
between Si x Ti and S 2 x T 2 - 

For a poset S and a subset I C S, an element t G S" is called an upper bound 
(resp. a lower bound) of X ii t > x (resp. t < x) for every a; G X . 

Definition 9.2. A poset S is called a lattice if for each pair of elements x,y, 
there exists the least lower bound x V y called supremum and the greatest lower 
bound x Ay called infimum. 

A lattice S is complete if every subset H C S (not necessarily finite) has 
supremum and infimum in S, and it is conditionally complete if supremum exists 
for each subset bounded from above, and infimum exists for each subset bounded 
from below. Note that the last two statements are equivalent so that, formally, 
only one of them is needed. 

Example 9.3. Consider the set of natural numbers N ordered in such a way 
that ni ^ 712 if and only if m divides n 2 . Then it can be verified that m V n 
is the least upper bound of m and n while m A n is their greatest common 
divisor. Other related examples are the lattice of subsets ordered by inclusion 
where A = (~l and V = U, or the lattice of convex sets ordered by inclusion where 
A = n but V is the convex hull of the arguments. However we emphasize that 
such examples are not important to us here. 

Example 9.4. Any linearly ordered set S, i.e. such that for each x,y G S there 
is x ^ y or y ^ x, is a lattice where both x V y G {x, y} and x A y G {x, y}. 
For example we can take the real line R or any subset of the real line, e.g., 
an interval [a,b] or a set of nonnegative numbers R + . However, R and R + 
are only conditionally complete, and for their completion we can consider R := 
R U {-00} U {00} and R + := R+ U {00}. 

From this one can construct slightly more complicated examples, which will 
not be linearly ordered, e.g., by means of Cartesian products. The following 
lattice-theoretic definition is of particular importance to us. 
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Definition 9.5. Let S be a set partially ordered by a relation ^. A (closed) 
interval in S is a subset of the form x := [x, x] = {t £ S: x ^ i ^ x}, where 
x, x e S and x ^ x. The elements x and x are called lower bound and upper 
bound of x respectively. 

Intervals can be viewed as pairs of lower and upper bounds. The set of such 
pairs is denoted by I(S) and called interval extension of S. It is evident that 
I(S) C S x S (but I(S) ^ S x S) and that I(S) is a poset. It can be shown 
that I(S x T) = I(S) x J(T) for posets S and T. 

Consider an algorithm A on posets. It takes input data (x\, . . . ,x n ) and 
generates output (y±, . . . , y m ). Here Xi £ Si for i — 1, ...,n and j/j 6 7} 
for j = 1, . . . ,m, where Si and Tj are posets. Hence, this algorithm induces 
a mapping 0^ : S —> T where S :— Si x • • • x S n and T := T\ x • • ■ X T m . 
We call an algorithm positive or nondecreasing if the corresponding mapping 
Si X • • • X S n — > Tx X • • • X T m is nondecreasing. 

Proposition 9.6. If an algorithm A is positive, then applied to the lower and 
upper bounds of an interval of S%X ■ ■ ■ X S n where Si , . . . , S„ are posets, it yields 
exact interval bounds on the application of A to the whole interval. 

Proof. These properties follow immediately from the positivity of A. □ 

It is clear that such algorithm and its interval extension have the same 
complexity. Using Proposition 19. 61 the interval extension 1(A) can be defined 
on the bounds of intervals only, and the result of 1(A) is a mapping from 
/(Si) x • • • x I(S n ) to /(Ti) x • • • x I(T m ). The complexity of 1(A) only doubles 
the complexity of A. 

We proceed with the following observations. 

Proposition 9.7. Cartesian product of positive algorithms is a positive algo- 
rithmm. 

Proposition 9.8. Let S be a complete or conditionally complete poset. Opera- 
tions (x,y) H> x V y and (x,y) H> x A y are positive mappings SxS->S. 

Definition 9.9. A mapping / : S — > T will be called interval regular or, 
briefly, i-regular if the following condition hold: for any interval x C S there 
exists yCT such that /(x) C y, y g /(x) and y 6 /(x). 

For example, consider R + := R + U {+00} and unary operation x 1— > x~ 
defined as x~ := x^ 1 on finite numbers, 0~ = +00 and oo~ := 0. This map- 
ping is not positive, but z-regular. Observe that an interval [a, b] is mapped to 
[b-\a-i]. 

Any algorithm A induces a mapping hence we can define i-regular algo- 
rithms. 

Definition 9.10. An algorithm A is called i-regular if f a is i-regular. 

Definition 9.11. An algorithm A is called ci-regularii fj± is interval regular and 
this algorithm and its interval extension have the same (asymptotic) complexity. 



20 



For example, the inversion operations in semifields and completed semifields 
are ci-regular but not positive. The following simple and obvious proposition 
illustrates a typical application of the notion of ci-regular algorithms. 

Proposition 9.12. Any composition or Cartesian product of ci-regular algo- 
rithms is a ci-regular algorithm. 

Corollary 9.13. Any positive algorithm is ci-regular. The complexity of its 
interval version is the same. 



9.3. Interval analysis over a fixed basic semiring 

Fix a basic positive semiring K , e.g., a basic numerical semiring in the sense 
of Subsection 7.1. For the sake of simplicity, we suppose that K is complete. 

So input output data run Cartesian products of several copies of K. Note 
that Cartesian product of positive semirings Si x 52 is a positive semiring, with 
respect to the Cartesian product of orders in S% and S%. 

Definition 9.14. An algorithm A is called elementary if it can be realized as 
a composition or Cartesian product of a finite number of basic operations from 
K, see Subsection 7.4 above. 

Algorithms which are not elementary might use if-else constructions, how- 
ever, they may be still ci-regular. 

From the results presented in Subsection 9.2 we can easily deduce the fol- 
lowing result. 

Theorem 9.15. The basic unary operation x 1— > x* induces a nondecreasing 
mapping S — » S , and the basic binary operations (£>, 0, as well as supremum 
and infimum, induce nondecreasing mappings S x 5* — > S . If an elementary 
algorithm V uses only these basic operations over a positive semiring, then it is 
positive. If K is a completed semifield, then the unary inversion operation is ci- 
regular but not positive. Every composition of elementary (and other ci-regular 
algorithms) is a ci-regular algorithm. 

Example 9.16. There exists an algorithm for computing A* which can be 
represented as a composition of positive operations with respect to K, see Sub- 
section 7.2 above. Hence this algorithm is positive. Moreover each standard 



algorithm solving the stationary discrete Bellman equation (see, e.g. [41|) is 
positive and ci-regular. 

Some standard algorithms for solving the matrix equations Ax = b, Ax = 
Bx,Ax = By, and other problems of tropical/idempotent linear algebra use 
the so-called (binary) pseudodivision operation, see, e.g. @,|2(J. In principle, 
thi s op eration is not elementary or ci-regular. However, using results presented 
in [20( , it is possible to reduce this operation to ci-regular algorithms for all the 
basic numerical semirings. So the corresponding algorithms of linear algebra 
are ci-regular. 
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